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We study the role of a possible nonet of light scalar mesons in the still interesting n —> 3n decay 
process, with the primary motivation of learning more about the scalars themselves. The framework 
is a conventional non-linear chiral Lagrangian of pseudoscalars and vectors, extended to include the 
scalars. The parameters involving the scalars were previously obtained to fit the s-wave tttv and 
7rK scatterings in the region up to about 1 GeV as well as the strong decay rf — > r/mr. At first, 
one might expect a large enhancement from diagrams including a light a(560). However there is 
an amusing cancellation mechanism which prevents this from occurring. In the simplest model 
there is an enhancement of about 13 per cent in the 77 — > 37r decay rate due to the scalars. In 
a more complicated model which includes derivative type symmetry breakers, the cancellation is 
modified and the scalars contribute about 30 percent of the total decay rate (although the total is 
not significantly changed). The vectors do not contribute much. Our model produces a reasonable 
estimate for the related ao(980) — /o(980) mixing strength, which has been a topic of current debate. 
Promising directions for future work along the present line are suggested. 



PACS numbers: 13.75.Lb, 11.15.Pg, 11.80.Et, 12.39.Fo 



I. INTRODUCTION 



There has been a revival of interest recently [|l] in the possible existence of a broad scalar meson (sigma) with a mass 
in the 560 MeV region and its corresponding nonet partners. A large number of workers p[-p8| have found evidence 
for the sigma in models of 7T7t scattering even though it is partially obscured by background. Generally this state is 
considered to be of exotic nature (more complicated than qq) and hence an important clue to an understanding of QCD 
in its low energy non-perturbative regime. Similarly, analyses of 7T7r , 7rK and tttj scattering have provided evidence 
for the existence of the remaining members of a possible light scalar nonet: the k, the ao(980) and the /o(980). In 
fact, the latter two states have been well established experimentally for some time. Of course, the treatment of such 
strongly interacting processes is inevitably model dependent and there are a number of different opinions as to the 
correct approach Thus it is of great interest to see whether treatments of the role of scalars in other processes 
using the same models employed in the scattering processes above give consistency with experiment. 

From this point of view we will study the role of possible light scalars in the interesting 77 — > 37r decay. Typically 
this process has been treated by chiral perturbation theory 29 , in which the possible effects of scalars have been 
amalgamated into effective contact interactions among the pseudoscalars. This is probably the most effective way to 
study the r\ — > 37r decay. However, our goal here is to learn more about the scalars so it is natural to keep them rather 
than integrating them out. Also there is a possibility that a light scalar [like the <r(560)] might give an enhancement 
due to closeness of its propagator to the pole [see for instance Feynman diagrams like (a) and (b) of FigJ|]. Another 
reason for including light scalars explicitly is to become more familiar with the isospin violating an(980) — /o(980) 
transition which should play a role in the 77 — > 37r decay and has also recently been postulated [g0| to provide an 
explanation for observations of anomalously strong ao(980) central production and the large T(cj> — ► /o7)/T(0 — ► ao7) 
ratio. It is important to know whether the value consistent with the eta decay determination is consistent with these 
proposed new effects. Doubts about whether an unreasonably large value was assumed in pfj were expressed in [try. 
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These doubts were confirmed |32j using the work of the present paper. Still another reason for the interest in the 
effects of the scalars in 77 — * 3-7T is to provide an orientation for the discussion of the apparently puzzling rf — > 3n decays 
in which light scalar mesons can be reasonably expected to have very large effects. We will give only a preliminary 
discussion of this process here. 

In section II we give a brief historical outline of treatments of rj — > 3n decay based on chiral symmetry. A number 
of well known ambiguities in the analysis are briefly described. 

Our calculation is based on the tree level treatment of a chiral Lagrangian containing pseudoscalars, vectors and a 
postulated nonet of light scalars. Since the calculation is somewhat complicated, it seems to us helpful to present the 
results in a series of steps. First, in section III we give the results of using a Lagrangian containing only pseudoscalars 
with minimal symmetry breaking terms. 

To this Lagrangian we add, in section IV, the scalar mesons. It will be seen that the individual scalar diagrams 
are quite large but there is a lot of cancellation so that the net effect is not at all dominant. However the scalars 
do, as desired, increase the predicted decay rate in a noticeable way. Next, the effect of adding some derivative type 
symmetry breakers for the pseudoscalars is described in section V. This doesn't much change the overall rate but 
modifies the somewhat delicate cancellations so that the scalars end up making a larger percentage contribution than 
before. In low energy calculations of this sort one always may expect some contributions from the vector mesons. 
This is discussed in section VI where it is shown that, although there is a new type of diagram the vectors do not 
produce a big change in the previous results. 

Section VII contains a discussion of the results and directions for further work. For the convenience of readers, 
material describing the chiral Lagrangian used is brought together in Appendix A. Similarly the detailed expression 
for the decay amplitude is given in Appendix B. 



II. HISTORICAL BACKGROUND ON THE 77 -> 3tt DECAY 



The study of 77 — > 3n has turned out to be surprisingly complicated and correspondingly important for understanding 
the non-perturbative (low energy) structure of QCD. Chiral dynamics in various forms has been the basic tool. Since 
the process violates G-parity it was initially assumed to be of electromagnetic nature, mediated by an effective photon 
exchange operator proportional to the product of two electromagnetic currents. The old "current algebra" approach 
had previously predicted the Kl —> 7r + 7r~7r° spectrum shape |33| to be 

1-^, (1) 

TO 

where to is the Kl mass and Eq the energy of the ir° in the Kl rest frame. This shape, which is in reasonable 
agreement with experiment, resulted from the vanishing commutator of the axial charge transforming like a tt + with 
the appropriate product of two weak currents. When Sutherland [j34| repeated this type of calculation for 77 — > 7r + 7r~7r° 
with the product of two electromagnetic currents he found that the decay amplitude was actually zero (to this leading 
order). Thus the 77 — > 3n decay did not seem to be mediated by a virtual photon emission and reabsorption. In fact, 
it was found p5[ that a quark scalar density operator with the AI = 1 property proportional to 

uu — dd (2) 

would give a non-zero result for the decay rate. A more detailed treatment p6| showed that the quark density operator 
gave the same spectrum for 77 — > 7r + 7r~7r as in Eq.(|l|) with to the 77 mass in this case. Such a result is in fairly good 
agreement with experiment. The scalar density interaction in Eq.(^) was recognized [p7| to be the fundamental 
up-down quark mass difference generated by the Higgs boson in the electroweak theory. 

However, the predicted rates of the 77 — > ir + ir~ir° and 77 — > 3n° modes (both the ratios and the absolute values) 
did not agree well with experiment at that time. Some years later, after more precise experiments, the ratio of the 
rates for 7r + 7r~7r° to Sir modes stabilized around the value expected from isospin invariance. On the other hand the 
absolute rate has only recently stabilized to a value notably larger than that predicted by theory. The theory behind 
the current algebra results could be economically presented in the framework of an effective chiral Lagrangian. For 
most low energy processes where the scheme could be expected to work, the tree level computation did produce results 
within 25 % or so of experiment. Thus the relatively poor prediction for 77 — > 37r at tree level is somewhat surprising. 



An improvement was obtained by Gasser and Leutwyler 1 29 who carried the computation of the chiral Lagrangian 



amplitude to one loop level. Since the non-linear chiral Lagrangian is non-renormalizable, this required the addition 
of new counterterms. Their finite parts were new parameters which could be mostly determined from other processes. 
They obtained the result T(r) —>■ 7r + 7r~7r ) = 160 ± 50 eV which may be compared with the present experimental 
value Q r (77 — > 7r+7r_7r °) cxpt = 267 ± 25 eV. The extra effects included involve both the implicitly integrated-out 
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heavier meson exchanges and partial unitarization to one loop order. One might expect a two loop calculation in 
the chiral perturbation scheme to be valuable but this may involve too many unknown parameters at the present 
stage. A dispersion approach using the Gasser-Leutwyler result as a subtraction gave an improved estimate |||] 
r (?y — > 7T + 7r~7r°) = 209 ± 20 eV, which still seems too small. 

A possible source of ambiguity arises from the determination of the coefficient of the driving scalar density interaction 
in Eq.(|^). This is determined from the K° — K + mass difference, which in turn has two components 



m 2 (/f°) - m\K+) = [m\K a ) m 2 (K+)] + [m 2 (K°) - m 2 (K + )] 7 , (3) 

corresponding to the quark mass differences and the virtual photon emission and reabsorption diagrams respectively. 
The latter is given in the chiral limit by m 2 (7r°) — m 2 (7r + ) according to Dashen's theorem |4(| and the reasonable 
assumption that the photon contribution saturates the pion mass difference. A number of authors f4l|| have argued 
that there are important corrections to Dashen's theorem which have the effect of boosting the 77 — s- 37r decay rate. 

If one questions Dashen's theorem it is natural to also question Sutherland's result, which deals with the direct 
electromagnetic contribution to 77 — ► 'in . An investigation of this point yielded Jl2"| the estimate that there was only 
about a 2% change arising from this, although it decreased rather than raised the rate. 

Still another point which may repay further investigation concerns the possible subtleties arising from 77 — 77' mixing. 
An understanding of the 77' — > 37r process, for example, might clarify this point. This process has been treated by 
some authors [Q , j44|] in the literature but has received only a fraction of the attention given to 77 — > 3-7T . 

In the present paper we will focus on learning more about the putative nonet of light scalar mesons by studying 
their contribution to 77 — > 37r. 



III. CHIRAL SYMMETRY RESULTS TO LOWEST ORDER 



For comparison, we first present the well-known results when only the terms present in the lowest order chiral 
Lagrangian of pseudoscalars are kept. 



LO . Elv (Vre ^ + s'Tr [M (U + *)] + ih" (|§L) 



(4) 



where the last term [see Eq.(A14) and comments there] supplies mass to the SU(3) singlet state and M. is defined in 
( All ). Fitting £ io to the experimental masses determines 8' — F^m^/8. 

The 77 — > 7T + 7r~7r° amplitude receives, in this approximation, contributions from diagrams (a), (b) and (c) of Fig]| 
which are given in Eq.(Bl) (with the non leading corrections deleted). To a reasonable approximation which displays 
the key dependences these sum up to the lowest order result for the 77 — > 7r + 7T~7r° amplitude 



Mo+-(E 1 ,E 2 ,E 3 ) 



(•5) 



Here Ei is the ir° energy in the 77 restframe and y is the dimensionless parameter in Eq.( All ) which measures the 
isospin violation in the quark mass matrix. Assuming Dashen's theorem, Eq.(H) yields 



Vy = Fl(m\ a - m 2 K+ ~ m 2 + m 2 + ), 



(6) 



whi ch al lows us to solve for y. Furthermore 9 P is the "nonstrange-strange" pseudoscalar mixing angle defined in 
Eq.( A15| ); it is generally taken to be about 37°. It is related to the "octet-singlet" angle, 9 by 



cos(L 



cos# — -\/2sin6> 
V3 " 



(7) 



Then Eq.(|5|) agrees with Eq. (1.14) of |29 except that they neglected 77 — 77' mixing by replacing cos6> p — > l/\/3 in 
what was denoted the current algebra formula. The matrix element for 77 — > 37r° is given in general by 



Mono — Mq + -(Ei, E 2 , E 3 ) + M ^ — (E 2 , Ei, E3) + Mq-| — (E 3 ,E 2 ,Ei). 



(8) 
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The widths are then, defining Ton = T(r] — > it + it 7r°) and Tooo = r(?7 — > 37r° 

1 

647r 3 m 
1 

3847r 3 m 



' 1 1 



To+- = g ,_;„_ / dE 1 dE 2 \M Q +-\ 2 , 



Fooo - / d^d^lMoool 2 . (9) 



Using £lo: with parameters determined as described above, we get the tree-level results [from the first terms in 
each of M^ act in Eqs.@]: 

r 0+ _ = 106 eV, 

^21 = 1.40. (10) 

These may be compared with the experimental results ]38| ] 

(r 0+ -) oxpt = 267 ±25 eV, 

(^521) = 1.40±0.01, (11) 
r o+- oxpt 

which demonstrate the disagreement with experiment for the overall rates in the simplest model. However the width 
ratio has about the correct magnitude. The related energy spectrum is also about the correct magnitude. The squared 
matrix element is usually described by quantities a, b and c defined from 

|M 0+ _| 2 oc {1 + aY + bY 2 + cX 2 ...), (12) 

with X = — — (E 2 — E 3 ) and Y = ^ — (E\ — m w ) — 1. In the present paper we shall not take into account the 

(not completely negligible) kinematic 7r° — tt + mass difference. Sec [^9| for a discussion of this point. The predictions 
from this simple model, a « — 1 and b « 0.25 are similar to the experimental results jl5) a oxp = —1.19 ± 0.07 and 
fo cxp = 0.19 ± 0.11 with c=0. 

It is of some interest to also give the predictions for the rj' — > 37r decay process at tree level using the simple 
Lagrangian Eq.ffl). It just is necessary (see Appendix B) to replace cos6 p by sin# p and m, q by m r y in Eq.(||) to get 
for the rl — ► 7r°7r~ 7T _ matrix element: 

M^Ek^Es) » ^sin0 p (l - ^). (13) 

r — Win' 



T' Q+ _ - 497eV, 



This leads to the predictions for the rf modes: 

T' 000 = 562eV. (14) 

The experimental results are given as |38| ] 

(rJ, + _)ex P <10 4 eV, 

(rooo)exp =315±56eV. (15) 

Only the 37r° mode has really been measured; its width is smaller than predicted in the simple model just presented. 
One would, of course, expect better agreement for the low energy process rj — > Stt for which chiral perturbation theory 
should be more clearly reliable. In the present paper we shall just make a few remarks on this more complicated 
process. 

It may also be worthwhile to give a rough estimate of the corrections to the rates corresponding to violations of 
Dashen's Theorem mentioned earlier. If we parameterize the electromagnetic contribution to the K + — K° mass 
difference as 

(m 2 K o - m 2 K +) y = /(m*o - m^+), (16) 

where / = 1 corresponds to Dashen's Theorem, we would find by using Eq.Q that the rj, rf — > 37r rates predicted for 
Clo should be multiplied by 

(m\ a - m 2 K+ ) - f{m 2 a - m£+ ) " 



(m% -m 2 K+ ) - (mlo 
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For / « 2, which was actually found many years ago |46|| the correction factor is about 1.54 and would give Tq^ « 163 

cV. This corresponds to the overall factor, y taking the value —0.33 while the Dashen's theorem value used to obtain 
Eqs.© and @ is y = -0.277. 




(c) 

FIG. 1: Feynman diagrams representing the pseudoscalar meson contribution to the decay rj — > 7r + 7r - 7r' 




(d) 



FIG. 2: Feynman diagrams representing the scalar meson contributions to the decay r\ — > tt + tt n 



IV. INCLUDING LIGHT SCALAR MESON INTERACTIONS 

Now we will study what effects the inclusion of a nonet of light scalar mesons will have on the 77 — > 37r calculation. 
We designate the scalar nonet by the 3x3 matrix N% whose interactions are also listed in Appendix A. N is assumed 
to contain the well-established /o(980) isoscalar and the ao(980) isovector as well as the cx(560) and the strange k(900). 
Of these only the k(900) will not contribute to 77 — > 3tt at tree level. The quark structure of such a nonet has been 
the subject of much discussion p|-p8[. If this were an ideal nonet like p — lu — K* — <f> one would expect the roughly 
degenerate eto(980) and /o(980) to be lowest rather than highest in mass. Actually the masses are better understood 
intuitively j47j] if N% is an "ideal dual nonet" constructed as Q a Q b with Q a ~ £ a bcq b q c ', qa being the ordinary quark. 
Then the observed inverted mass ordering is easily seen to follow just from counting the number of strange quarks in 
Ng. It is important to note that the form of the couplings of N% to the particles of the nondinear chiral Lagrangian 
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(a) (b) 

FIG. 3: Feynman diagrams representing the p meson contributions to the decay 77 — > 7r + 7r~7r° . 

being used depend only on the flavor transformation properties of N%. This does not distinguish different quark 
substructures. What is sensitive to the quark substructure is the scalar mixing angle, # s , defined from 




(18) 



Small values of 9 S would typify a dual ideal nonet while \9 S \ about ^ would typify a conventional nonet. Fitting 
the 7T7T and irK scattering amplitudes, including the effects of these scalar resonances, selects [|l6| the small value 
9 S = -20.3°. 

The scalar nonet mass terms in Appendix A are specified by the four parameters (a, 0, c, d). The needed chiral 
invariant scalar -pseudoscalar-pseudoscalar S<p(p-type couplings are specified by the parameters (A, B, C, D) and the 
mixing angle 9 S . These were all determined from fitting to tttt scattering, 7rK scattering and the strong decay 77' — ► rjirir. 

Actually there is reason to believe [Q that the scalars may be best understood as mixtures of a lighter dual nonet 
and a heavier ordinary nonet. From that point of view, which will be explored more fully in the future, the present 
single nonet, N£ should be regarded as an approximation to the situation where the heavier (after mixing) particles 
have been integrated out. 

The Feynman diagrams for the scalar contributions to rj — > 7r + 7r~7r° are shown in Fig.^J. Notice that the diagram 
in (d) involves new ao — a and ao — /o isospin violating transitions rather than the tt° — r\ and tt° — 77' transitions 
which play an important role in the other diagrams. Their strengths A aa and A a f (see Appendix B) were determined 
simply by including the effects of isospin violation contained in the spurion matrix M. in the b and d scalar mass 
terms. Therefore, this does not introduce any new parameters. Actually, the possibility of such contributions was 
suggested a long time ago |43|| as a possible solution of the rj — > 3n width problem. Recently a relatively large ao — fo 
mixing has been suggested ]30( | as a way of understanding both anomalously large ao central production and the large 
T((j) — > /o7)/r(</> — > ao7) ratio. However criticisms of this explanation have also been presented Q. Clearly it may 
be useful for studies of processes other than 77 — > 37r to give the coefficients of the scalar isospin violating two point 
Lagrangian, 

£ = A aa a°<r + A af a° f , (19) 



determined consistently with the 77 — > 37r calculation. Using the parameters from Eq.(A21) in Eq.(B4) we find 
Aaa = 0.0I70y GeV 2 , A af = 0.0234y GeV 2 , where y is the quark mass ratio ™ u +™ d . Notice that y (which is 
negative) is an overall factor for the 77 — > 3n amplitude in the present model. 

We would like to discuss the effects of the scalars when added to the more realistic Lagrangian presented in Appendix 
A which contains both pseudoscalars and vectors. This Lagrangian contains additional symmetry breaking terms (a p 
and A') to account for the ratio of pseudoscalar decay constants, F^fF^ being different from unity as well as a number 
of terms describing the properties of the vector mesons. Of course, the vector mesons play an important role in low 
energy processes. Since there are many terms it seems useful to add these new features one at a time. Thus in the 
present section we will consider just the minimal pseudoscalar Lagrangian, C^o [Eq.@], to be present in addition to 
the scalars. Furthermore, it is instructive to look at the contributions to the amplitude from different diagrams in 
order to see how they combine to give the predicted total 77 — ► 7r + 7r~7r° width. In section III we reviewed the leading 
order calculation of the 77 — > 37r amplitude which gives the result r(?7 — > 7r + 7r~7r°) = 106 eV in the Dashen's theorem 
limit. In Fig.|| we show how the individual contributions of the diagrams in Figjl combine to give the leading order 
amplitude. The magnitude of the 77 — 7r° and 77' — 7r° transition coefficiencts [Eq.(B4) with a p = A' = 0, so independent 
of which state is on- shell] are (in GeV 2 ): 

C nri w 0.0042, « 0.0031. (20) 
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FIG. 4: Plot of different contributions to the leading order r\ — > 7r 7r - 7r amplitude as a function of the energy of the neutral 
pion. On the left, the solid line corresponds to the direct four-point isospin-violating r]ir + -ir~Tv° vertex of Fig.hla. The energy- 
dependent dotted line is due to an r/iv transition followed by a four- pion contact vertex (Figjljb). The dashed line is the small 
contribution due to Figjljc. (with both rj and rj included in the final line). On the right we show the total leading order 
amplitude, which is the sum of Figs.nla-lllc. 



Three of the twelve scalar diagrams in Figj2| can be seen to be larger (by at least an order of magnitude) than all 
the rest. These are the three diagrams involving the lightest of the scalar mesons, a, and are contained in Figs.|]a, ^3 
and |d. For the case of Fig.||b the graph with an r\ — ir° transition dominates that with an rj' transition because the 
latter is suppressed by the (square of) the rf mass in the denominator of the propagator and also the smallness of the 
associated coupling constants/transition coefficients. In Fig.^ we present the rj — > 7r + 7r~7r° amplitudes arising from the 
three diagrams just mentioned and notice that they cancel almost completely. In particular the a exchange diagrams 
in Figs.||a and gb have opposite signs, as expected - their structure is roughly similar except the propagators have a 
relative minus sign. We note also that the new isospin violating diagram, involving an ao — u transition, turns out not 
to lead to dramatically larger contributions than the other diagrams. The cancellation between different diagrams 
involving the sigma means that the total scalar contribution to the rj — > 7r + 7r~7r° width is smaller than might be 
expected and in fact arises mainly from the a,Q exchanges in Fig.^c. Comparing Fig.^Jb with Fig.^b shows that the net 
scalar contribution does enhance the overall 77 — * 7r + 7r~7r° rate. Specifically, including the scalar contributions with 

the pseudoscalar Lagrangian Ceo has increased Ton by 16 per cent, from 106 eV to 124 eV. The ratio Tooo/ron 

is essentially unchanged. 

Actually, the calculations above have neglected the finite widths of the a, fo and ao particles. We take these into 
account by making the replacements in the corresponding propagators [see Eq.(B2)]: 



: > : j (21) 

m x + I 2 m x + I 2 ~ i m x^x ' 

where X stands for a, fo or ao. The Tx are given in Appendix A. These replacements modify the result to Tq^ = 120 

eV, a 13 per cent increase relative to the leading order result. The width effect is mainly due to the a propagator. 

We may note that the improvement due to the scalars is consistent with the lower values of the prediction, 160 ± 50 
eV obtained from the next order of chiral perturbation theory in rcf . [p9[ . The numerical amount of suppression of 
the scalar contribution to t he dec ay rate from cancellation of Figs, (a) and (b) of Fig.|| is due to the fitted values of 
7 CT7rTr and 7 CTr;l) given in Eq.(A22). If we wanted to raise the predicted rate to about 150 eV (still keeping Dashen's 



theorem in the evaluation of y) it would be necessary to raise ^a-q-q to about 10. 



V. EFFECTS OF HIGHER ORDER PSEUDOSCALAR SYMMETRY BREAKERS 



So far we have worked only with the leading order chiral Lagrangian of pseudoscalars and sc alars and obtained (to 
linear order in y = — " r l d +™ u ) isospin- breaking amplitudes proportional to S',b and d in Eq. ( |A13| ) of Appendix A. 
In order to better fit the propert ies of the pseudoscalar mesons we consider, as mentioned above, the higher-order 
symmetry breaking terms in Eq.(Al3) with coefficients a p and A'. The numerical values of these parameters are 
obtained in section 4 of Appendix A, based on an overall fitting of pseudoscalar meson properties. 
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FIG. 5: Scalar meson contributions to the 77 — > n ir~n amplitude as a function of the energy of the n . On the left we show 
the individually largest contributions - the solid line corresponds to Fig.|^, with a exchange, the dotted line corresponds to 
Fig.^b with a exchange and the 77-71-° transition while the dashed line corresponds to Fig.^d involving the new isospin violating 
ao — <y transition. On the right we show the total amplitude due to the scalar mesons alone. The solid line is due to all of the 
diagrams in Figj^ (for the sample value E2 = m,r) and the dashed line is the sum of the three largest amplitudes plotted on 
the left and discussed in the text. 



Next we examine the effects of these two new symmetry breaking terms on our previous calculation. It will be seen 
that these effects include an interesting redistribution of the contributions from the scalar and pseudoscalar diagrams 
to the total amplitude. First, the contact diagram Fig.[l]a will receive corrections due to the a p and A' terms [see 
Eq.(Bl)]. This results in some energy dependence since a p gives a four-point derivative coupling. More importantly, 



the 77 — 7T° and 7/ — ir° transition coefficients relevant for 77 — ^ 37r now depend on which particle is on-shell and are 



numerically (in GeV 



ryn 



-0.00583y, 

TTTj f 



^ '7T7] 

-0.01132/. 



-0.0151t/ 



(22) 

Since C% v is now considerably suppressed in magnitude, the Feynman amplitude for Fig.|l|b is now suppressed, while 
CJ I( and the amplitude for Fig.^3 remain about the same. These results, due to only pseudoscalars, are summarized 
in Fig.^J which may be compared with FigJ^. The net result is that the total 77 — > 7r + 7r~7r° amplitude (shown in the 
second of FigJ|) due to pseudoscalar mesons is reduced compared with the leading order result. The pseudoscalars 
themselves now give To-i = 81 eV rather than 106 eV, as in section IV. 



-0.05 - 



0.177 
E1 GeV 



0.157 0.177 
E1 GeV 



FIG. 6: Left: Contributions due to pseudoscalar mesons alone, taking into account higher order symmetry breaking effects 
encoded in a p and A'. Compare with the first of Figs.^J Right: Total amplitude due to pseudoscalar mesons alone, taking into 
account a p and A' terms (compare with second of Figs.^]). 

However, for the diagrams involving scalar mesons the effect of higher order symmetry breaking is even more 
important. As we noted above, the scalar meson contribution to 77 — > 7r + 7r~7r° was rather small as the main diagrams 
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tended to cancel. When we include the a p and A' corrections to the -qir transition this cancellation will not be so 
complete. Specifically, comparing Eqs.(pO|) and ( |22| ) we see that the magnitude of the amplitude for Fig.||a, where the 
77 — 7T° transition occurs with an on-shell 77, will be reduced by a factor of approximately four, while that of Fig.Hb, 
where the rj — ir transition has an on-shell pion, will remain about the same relative to our result in Section IV. Fig^p 
will be unchanged. There will now be a non-negligible contribution from the scalar mesons. It will be more negative 
in sign (the contribution from FigJ^a is positive, but now smaller in magnitude) and will add "constructively" to the 
pseudoscalar diagrams in Figjl] and so increase our prediction for T(rj — > tt°it + tt~). This is shown in Fig. ^. Adding 
all of the pseudoscalar and scalar diagrams in Figjl] and FigJ| with the inclusion of the symmetry breaking effects due 
to a p and A' we get: r(ry — > 7r°7r + 7r~) = 119.6 eV. This is essentially the same as the result in section IV but now a 
larger portion is due to the scalar meson diagrams. Since this is the case the damping effect of the sigma width will 
be more prominent; in fact it reduces the predicted rate to 103.6 eV in the Dashen's theorem limit. 




E1 GeV E1 GeV 



FIG. 7: Left: Contributions due to the three largest scalar meson diagrams, taking into account higher order symmetry breaking 
effects encoded in a p and A'. Compare with the first of Figs^. Right: Total amplitude due to pseudoscalar and scalar mesons, 
taking into account a v and A' terms (compare with second of Figs.^]) . 



It may be of interest to see how the detailed pattern just described depends on the precise value of the quark mass 
ratio x and, as explained in section 4 of appendix A, correspondingly on the crucial isospin violating quark mass ratio 
y. This is shown for the predicted value of Tq-i in Table |[ 



x, y ps. only ps. + sc. 


(zero scalar widths) ps.+sc. (non-zero scalar widths) 


20.5, -0.202 63.7 eV 95.9 eV 


82.2 eV 


23, -0.241 70.7 eV 106.0 eV 


91.4 eV 


25.1, -0.277 80.2 eV 119.6 eV 


103.6 eV 



TABLE I: r(ry — > -k + tv~-k { ) for different values of x and y defined after Eq.(All). The second column applies to the case 
of only pseudoscalars present while the third includes scalars too. The effect of taking non zero scalar widths into account is 
shown in the last column. Dashen's theorem is assumed in order to extract y. 



VI. INCLUDING VECTORS IN THE CALCULATION 



It is well known that the inclusion of vector mesons is important for a realistic discussion of low energy chiral 
dynamics. For example, in the chiral pertubation scheme, most of the finite pieces of the counterterms can be 
explained by integrating out various vector contributions f49|| . In our present approach, of course, we are keeping the 
resonances, rather than integrating them out, in order to learn more about the scalars. 

First, the vector mesons contribute to the 77 — * tt + tt^tt° amplitude corresponding to Fig.^ja, which is just a correction 
to 7T7r scattering. Its value given by the amplitude, in ( |B3 j ) is easily seen to be comparatively large. However 
the fourth term of the U(3)l x U(3)r invariant Lagrangian Eq. QAlO] ) gives, in addition to the pirir vertex a four pion 
contact term [included in M® ontact in (Bl)]. Actually this contact term cancels most of the contribution from the 
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p-exchange diagram in Fig.|3|a. This is well known from chiral treatments of irir scattering: when the p is added to 
the Lagrangian, chiral symmetry requires a contact term which cancels most of the p contribution near threshold, 
thereby maintaining the current algebra threshold result. 

However, the situation is actually a bit more complicated since the process of obtaining an adequate fit to the 
properties of both the vectors and pseudoscalars fl50| requires a number of additional symmetry breaking terms 
shown in Eq.( |Al3| ) of Appendix A. As well as the symmetry breaking terms we have already discussed involving the 
pseudoscalars alone, there are, in particular, two new terms, measured by the coefficients a + and a_. It turns out 
that their effects are very minor. They include an additional contribution to the 4-point isospin violating ry7r + 7r _ 7r° 
vertex due to the a_ term, corrections to the 4-pion vertex in Fig.[l]b due to both a + and a_ and an additional 
diagram, shown in FigJ^b, which contains a new G-parity (and isospin) violating piri] vertex [the amplitude for this 
is given as M h p in ( |B3[ )] . Note that there there exists a/)°-w mixing transition, which is the analog of the tt° — r\ and 
a o — fo mixing transitions, but it does not contribute to r\ — > 37r at tree level. 

The decay widths with inclusion of vectors are tabulated in Table || for the same values of x,y used in the last 
section. In this table the neutral modes are also included. Furthermore, the effect of both the scalar and (actually 
negligible) vector widths are included too. We see that, as expected from our discussion above, the vectors do not 
change the overall predictions compared to the last column of Table | very much but they do give a little enhancement. 
This is also clear by comparing the pseudoscalars + vectors column of table || with the pseudoscalars only column of 
Table B. It is seen again that the scalars make a non negligible contribution to the total amplitude. 



x, y decay mode ps. +vec. 


ps.+sc.+vec. (no 


width) ps. +sc. + vec. (width included) 


20.5, -0.202 0+- 


64.4 eV 


96.6 eV 


82.8 eV 


000 


92.9 eV 


139.4 eV 


118.9 eV 


23, -0.241 0+- 


71.9 eV 


107.4 eV 


92.5 eV 


000 


101.9 eV 


152.9 eV 


131.1 eV 


25.1, -0.277 0+- 


82 eV 


121.7 eV 


105.4 eV 


000 


114.5 eV 


171.3 eV 


147.7 eV 



TABLE II: F(r/ — * 7r°7r + 7r - ) and T(r) — > 3ir°) for different values of x,y. In the third column pseudoscalars and vectors are 
both present. In the fourth and fifth column pseudoscalars, vectors and scalars all present (without and with the effect of the 
scalar meson widths). 



Finally it is interesting to display the energy spectrum parameters a, b and c defined in Eq.dl2) for the various 
models we have examined. These are given in Table III and are seen to be reasonable. The \ 2 measures the fit of our 
model to the spectrum shape assumed in Eq.([12|) and seems to be small. 





a 


b 


c 


x 2 








pseudoscalars (LO) 


-1.11 


0.31 





5.6 


X 


10" 


-5 


pseudoscalars 


-0.96 


0.23 





1.5 


X 


10" 


-4 


pseudoscalars+scalars 


-0.93 


0.22 


-0.01 


3.3 


X 


10" 


-4 


pseudoscalars+scalars+vectors 


-1.09 


0.26 


0.033 


5.8 


X 


10" 


-3 



TABLE III: Fits of the energy dependence of the normalized (charged) decay amplitude for y — > ty°-k + -k~ to the form |Mo+_ 2 = 
1 + aY + bY 2 + cX 2 . The first line corresponds to result at leading order with pseudoscalar mesons only. The second with 
inclusion of higher order symmetry breakers, the third when scalar mesons are added and the final line when vector mesons are 
included as well. 



VII. DISCUSSION 



We studied the role of a possible nonet of light scalar mesons in the still interesting 77 — > 3tt decay process. Our 
motivation was primarily to learn more about the scalars themselves. The framework is a conventional non-linear 
chiral Lagrangian of pseudoscalars and vectors, extended to include scalars (the Lagrangian is described in Appendix 
A). The parameters involving the scalars were previously obtained to fit the s-wave irir and irK scattering in the 
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region up to about 1 GeV as well as the strong decay rj — > ijtttt. An initial concern is whether the model as it 
stands, containing essentially no undetermined main parameters (up to possible uncertainties in the quark mass 
ratios x = 2m s / (m u + md) and y — (m u — md)/(m u + 1x14)), does not make the 77 — > 37r amplitude too large. 

In particular, the cr(560) exchange diagrams (a) and (b) of Fig.|] might lead to a great deal of enhancement due 
to the possibility of the <r(560)'s momentum being close to mass shell. However this turns out not to be the case. 
In our initial calculation where the scalars are added to the minimal model of pseudoscalars given in Eq.(|), the left 
part of Fig.|] shows that these two diagrams, though not individually small, tend to cancel each other. This partial 
cancellation occurs because the 77 — tt° transition leads to opposite signs when the 77 is on mass shell and when the ir° 
is on mass shell (In the first case we have a ir° propagator carrying the momentum squared of an on-shell 77 while in 
the second case, the reverse holds). In addition, the enhancement due to the sigma propagator is further suppressed 
by the inclusion of an imaginary piece, needed to satisfy unitarity in the scattering calculation. The net result is that 
the effect of including scalars in the minimal pseudoscalar Lagrangian, Eq.(Q) is to increase the width for 77 — * 7r + 7r~7r 
decay by about 13 per cent. This relatively small, due to cancellation, increase illustrates the difficulty of finding 
dramatic "smoking gun" evidence for the existence of a light sigma. In the scattering calculation a light sigma appears 
(see for example |13| ) obscured by a large background and does not have a simple Breit Wigner shape. 

It is amusing to note the effect of higher derivative terms in the Lagrangian of pseudoscalars (see section V for 
details). The higher derivative terms allow one to conveniently implement at tree level the fact that the ratio of 
the pseudoscalar decay constants Fxp/F^p is somewhat greater than unity. With these terms the important ir°r] 
transition vertex has a momentum dependent piece. Together with a needed modification in the parameter fitting 
(see section 4 of Appendix A) this reduces the contribution of the pseudoscalars to the 77 — ► 3ir decay width. However 
the modification of the ir° — 77 transition noticeably upsets the cancellation between the two sigma exchange diagrams 
in (a) and (b) of Fig.||. The net result is that, while the total prediction for the 77 — s> 37r decay width remains about 
the same, now about thirty percent of the value is contributed by the scalars. 

The vector meson contribution, discussed in section VI, actually does not change things much. This is because 
the p exchange diagrams for tttt s-wave scattering are essentially canceled at very low energies by an extra four pion 
contact term which automatically arises due to the chiral symmetric formulation. Experimentalists fit the Dalitz plot 
describing the 77 — ► 7r + 7r~7r spectrum to the form given in Eq.([l2|). A fit of this type to the predicted spectrum 
from the Lagrangian of pseudoscalars, scalars and vectors was seen to be close to the experimental one. The basic 
spectrum shape is already reasonable with the very simplest model discussed in section III. As both the theory and 
experiment get more precise, the importance of the spectrum shape toward a deeper understanding of the underlying 
physics increases. 

A particularly interesting scalar contribution to 77 — > 37r arises from the ao — a transition shown in (d) of Fig.|^ 
(The ao — /o transition contribution to 77 — > 3ir is suppressed due to the propagator of the heavier /o(980)). This 
is the analog of the important 7r° — 77 transition and, in a sense, is a new mechanism for 77 — ► 3ir (although it was 
investigated a long time a go |43| as a possible way to increase the 77 — > 3tt width). The formula in raw form for this 
transition is given in Eq.(B4). We evaluat ed it s strength from the knowledge of the isospin violating piece of the 



dimensionless quark mass matrix M. in Eq.(All), determined from the pseudoscalar sector and the coefficients: a, 6, c 



and d of the scalar meson mass terms [see Eqs.( |A10| ] and ( A13)) determined from the isospin conserving sector of the 
scalars. However as one can see from the left sides of Figs. I and @, the contribution to 77 — > 3-7T due to the a - a 
transition is not very large, although it has the right sign to boost the decay rate. 

The method just described also evaluates the strength of the ao(980) — /o(980) transition. For convenience this 
is given after Eq. ([l9|) , where the overall factor, y is displayed. This transition has been very much "in the news" 
recently as a proposed [Q explanation for the large observed T(4> — ► /o7)/r(</> — * 007) ratio and the anomalously 
strong ao central production. However criticisms of this explanation have been given ]3l| ], p2| ], pointing out that 
the ao — /o mixing expected from a transition strength like the one determined above is insufficient to give a large 
effect. Intuitively, because of the near degeneracy of the ao(980) and /o(980) as well as the similarity of their widths, 
one might expect the mixing to be very large. But the mixing amplitude is governed by a dimensionless factor 
iA a f /(m a r a ) [see for example Eq.(12) of ^]] which is suppressed by the scalar meson width, r a . 

In section II we discussed the current comparison between theory and experiment for the 77 — > 7r 7r + 7r~ width. The 



experimental width 1 38 is Tq^ = 267 ± 25 eV. This may be compared with the one loop chiral perturbation theory 

result |2^] of 160 ± 50 eV. More recent attempts ||^] to estimate final state interaction effect outside of the chiral 
perturbation theory approach have increased this somewhat to 209 ± 20 eV. It seems to us that the thirty per cent 
contribution of the scalars compared to the pseudoscalars we have found should probably not be considered on top of 
this latter figure. That is because a good portion of the increase due to scalars we have found may be considered as 
resulting from final state interactions. Many attempts to close the gap between theory and experiment have focused 
JllJ on a reanalysis of electromagnetic corrections to the K + — K° mass difference. This is argued to increase the 
quark mass ratio, y which is an overall factor for the 77 — s- 37r amplitude. 

From the standpoint of learning more about the properties of the scalar mesons it is clear that the 77' — * 3ir 
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decays represent a potentially important source of information. In this case there is sufficient energy available for 
the a (980) and / (980) propagators to be close enough to their mass shells to avoid suppressing the contributions of 
these resonances. On the other hand, the theoretical analysis is more difficult since large non-perturbative unitarity 
corrections are expected. In addition, other more massive particles may also contribute. The experimental information 
[see Eq.(fn^)] is more preliminary than in the 77 — j- 37t case. While a number with reasonably small errors has been 

presented for Tqqq, there is only a weak upper bound for and also no information on its Dalitz plot. In the 

model employed in the present pap per the rj' — > Sir amplitudes are simply obtained from the 77 — s- 37r amplitudes by 
the simple substitution given in Eq.(|B5|). Notice that this substitution rule would get modified if a more complicated 

)] is adopted. As shown in Eq.(|l4|) the prediction of the 
minimal model of only pseudoscalars is somewhat too high, but at least of the correct order of magnitude. Adding the 
scalars without any readjustment of parameters does not improve the prediction for Tqqq but makes it considerably 
larger (about 2300 eV). A similar large value was recently found in H]. Since the phase space is fairly large it is 
perhaps to be expected that large values are typically obtained. Presumably it is a sign for including more detailed 
unitarity corrections or other physical effects which result in cancellations. We are particularly hopeful that a careful 
study of mixing between a lower mass exotic scalar nonet and a more conventional higher mass scalar nonet (2^, ^8], ^] 
may solve this problem and perhaps contribute to an improved understanding of the i] — > 3ir decays also. 



r) — rf mixing scheme [e.g. the one mentioned after Eq.(A15 
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APPENDIX A: THE CHIRAL LAGRANGIAN 

For convenience we collect here needed terms from the pseudoscalar- vector chiral Lagrangian presented in |5(| and 
from the scalar addition presented in ]lq]. 



1. Transformation Properties 



These are constructed to mock up the symmetry properties of the fundamental quark Lagrangian, under which left 
and right projected light quark fields tranfsorm as 

Qlm -> U L , R q L , R , (Al) 

Ul and Ur being 3x3 constant unitary matrices. The pseudoscalar nonet <ft(x) is a 3 x 3 matrix which fits into the 
unitary chiral matrix 

U = cxp(^^) (A2) 

where F„ is the (bare) pion decay constant. Under a chiral transformation 

U -> U L UU R . (A3) 
It is convenient to define the 3x3 unitary matrix £ by U = £ 2 . Then £ transforms as 

£ -» U L Zrt(4>, x) = K{<j>, x)£U R , (A4) 

which implicitly defines the unitary matrix K . The intuitive significance of K is that the objects Kq behave like bare 
quarks surrounded by a pseudoscalar meson cloud, or "constituent quarks" . The objects 

VvPv = ktd li ?±?d li Z), (A5) 
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transform as 

u„ - Kv^ + iK8^Kl 
A putative scalar nonet matrix iV(a;) is taken to transform as 

N ^ KNK\ 

The vector meson nonet p M transforms as 

p„ -> jr^t + I^A-t 
9 

where we have included the dimensionless coupling constant, g. The "field-strength tensor" 

FfiAP) = ®nPv ~ ®vP» ~ W \Pm pA -> kf^kK 



(A6) 
(A7) 

(A8) 

(A9) 



2. t/(3)z, x U(S) R Invariant Terms 



These comprise the kinetic terms for the three multiplcts, mass terms for the scalars and vectors and appropriate 
interaction terms: 



Co 



^-TT(d^d^) - ±Tr(F^(p)F^(p)) 



- -Tr^NV^N)- ^Tr [(gp, 

- aTi(NN) - cTr{N)Tr{N) 



+ Fl 



Ae abc e Aef Ni{pX{ Vil ){ + BTr(N)Tr(p^) + CTi(Np^Tr(p^ + DTv(N)Tv(p^Tv(p^] , (A10) 



where V^N = d^N — iv^N+iNv^. These include the parameters ml, a, c, A, B, C and D. Note that the pseudoscalars 
are still massless at this level. Further note that for the interactions and mass terms of the scalars we do not restrict 
ourselves to a single trace. For qq mesons the single trace is suggested by the OZI rule while for an ideal dual nonet the 
A term is in fact expected to be dominant. We made a fit for ml, a, c, A, B, C and D assuming only SU(3) invariance. 



3. Symmetry Breaking Terms 



The fundamental QCD Lagrangian contains the quark mass term — ( m "+" ld ) qJ\Aq, with the dimensionless matrix 



M 



1 + y 
1-yO 
Ox 



where x 



and y 



rn u +m d 

It is convenient to define 



rn d —ra u 
m d +m u 



Then, the symmetry breaking Lagrangian is taken as 

C SB = S'Tr [M{U + U^)] + X' 2 Ti [MU^ MU^ +MUMU] 



2a_ 



M- 



2a+Tr 



+ 2 7 'Tr \M+F^(p)F^(p) 

- bTr(NNM) - dTr(N)Tr(NM). 



(All) 



(A12) 



(A13) 
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Only the parameters 5' , A' 2 , a p , b and d here contribute to the isospin violating vertices of interest in the present 
paper. The parameter 7' was included in the overall parameter fit obtained in fr50|| but its small effect on the isospin- 
conserving vertices will be neglected. 

In addition to the quark mass induced symmetry breaking terms there is an important term induced by instanton 
effects which breaks just the U(1) A piece of SU(3)l x SU(3),r x U(l)y x U(l)yi- It may be summarized as 

k , o ( detZJ A / \ 1 a \ 

£ ^ 576 ln (dWtj + --- (AU) 

where k is a constant essentially proportional to the squared mass of the rj' meson. The three dots stand for other 
terms which will be neglected here but are listed in Eq. (2.12) of ]5l| ]. Effectively this term gives an important 
contribution to the rf mass and an 77 — rj mixing angle defined by 



?y' / \ sin6» p cos6> p J I <j>l 



(A15) 



When the extra terms in Eq. A14 are included they will not only give rise to an additional isospin violating transition 
but will also modify the 77 — rj mixing transformation above to be the non-orthogonal one given in Eq. (4.9) of |l). 
We will not include these effects in the present paper, however. 

4. Numerical values of parameters used 

For the averaged pseudoscalar masses we used, 

m„ = 0.137 GeV, m K = 0.4957 GeV. (A16) 

In section III we gave the fitted parameters for the lowest order Lagrangian containing only pseudoscalars. This 
also yields the isospin conserving quark mass ratio x — 25.1 (assuming that /, defined in Eq. ( [l6| ) is unity). A 
refitting of these parameters is necessary when the a p /g 2 and A' 2 symmetry breaking terms are included. This can 
be conveniently done following the method used in preparing Table III of |5Q ]. There, a value of x is assumed and 
the four quantities F v (unrenormalized pion decay constant), <5', |A'| 2 and a p /g 2 are calculated in terms of the four 
physical quantities m^, ttik, F^p = 0.1307 GeV and Fk p = 0.1598 GeV, using: 



A 



12 



(l + x)F 2 p ml/m-F 2 p m 2 K /8 



1 

= F 2 p m 2 /8-4A' 2 , 
2 



{Fk p I Fyp) — 1 



5 2 F 2 2(1 + x)-A{F K p/F^ P Y 1 

F n 



± up 



a 



if = n^). (ait) 



(l + 4a p /(.g 2 F 2 ))V2 

In addition, the isospin violating quark mass ratio y is obtained from 

(m 2 K0 - m 2 K+ ) - /(m 2 - m 2 + ) = (4y/F 2 p )(-25' - 8(1 + x)\' 2 + m 2 K a p /g 2 ), (A18) 

for a particular value of /. To isolate the effects of the scalars we may choose an x such that, with the value / = 1 
corresponding to Dashen theorem, we recover the value y — —0.277 found in sections III and IV. That gives x = 25.1 
and 

Fk = 0.128 GeV, 5' = 0.0386 x 10" 3 GeV 4 , 
a p /g 2 = 0.176 x 10~ 3 GeV 2 , |A'| = 0.643 x 10~ 3 GeV 2 . (A19) 

The needed dependences on the quark mass ratio ,x of the parameters involving vector mesons (7' , a + , a_, m v 
and g) are given in Table 3 of Q; the additional point x — 25.1 used in Table || above was treated by interpolation. 
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The masses and widths of the scalars are taken to be (in MeV) 

to ct = 550, m K — 897, m aQ = 983.5, m/ = 980 

T ff = 370, T ao = 70.0, r /o = 64.6. (A20) 

Note that the values of Y a and T K are not "Breit-Wigner" widths but a re cho sen t o unitarize the irir and ttK scattering 
amplitudes. The masses above fix the parameters (in GeV 2 ) in (A1C) and (A.13) 



a = 0.492, b = -0.00834, c = -0.0160, d = -0.00557 (A21) 

and the mixing angle 8 S — —20.3°. The parameters A, B, C, D define all the trilinear Scfxj) coupling constants according 
to the formulas given in Appendix C of p6f . The needed coupling constants are (in GeV -1 ) 

7<T7r7r = 7.27, j avv = 3.90, j ar}r} > = 1.25, ja V 'n' = -3.82, 
~f flT ir = 1.47, 7 /w = 1.50, 7/ w / = -10.19, jf n ' V ' = 1-04, 

lann = -6.87, 7w = -8.02. (A22) 

APPENDIX B: DECAY AMPLITUDE 

The Feynman diagrams representing the T)(p) — > ir° (pi)ir + (p2)tt~ (P3) decay are shown in Figs. 1-3. The contact 
diagrams (la, lb and lc) receive contributions from the pseudoscalar and vector part of the Lagrangian 

a™ .l6yS'cos9 p .8ya p cos6 p .512y\' 2 cos6 p 

M contact = 1 + ' 1~2 T?4 (~ 3 P2-P3 + P-Pl + P-P2 + P-Ps) + « ^4 " 

or^ p r np 7rp 

h ( 3to 2 \ C r > 2 . 

M CO ntact = +» I 1 ~ 4 - 2j?2 J m 2 _ m 2 ^pT {~^P2-PZ + Pi-Pi ~ P-P3+P1-P2 ~ P-P2 + 2p-Pl) 
+ 8 -2 F 4 T( ~ Z\\- 2 P-Pl + 2 P2-P3+P3-P ~ P3-P\ + P3-P ~ P2-P\) 

■ 8a,, Cg, 

+ ^^(TO^-TO^) ( 5P ' Pl ~ ~ ^ 3 + - ^ + ^ 



256A ,2 C*IL 16<J' C2 



contact 




P2 -P3 + P-P3 + P-P2 ~ Pi -P3 - Pi -P2 ) 



(Bl) 



The scalar contributions (FigsJ^a, b, c, and d) are: 



,,„ 2C; 7g ff „ (p.pi)(p2.p 3 ) 

M scalar = f 2 , / ^ + \ a ^ fo) 

m% — m„ + (p — p\Y 

,W, . ft I ICl lan-nlvrm C^'TowTawM {P-Pl){P2-P3) , , , , 

M -"- = ^ m ?) -TO2 + TO 2 / - TO^ J to2 +(p - Pi ) 2 + 

M c = ■( C nr^ V C^,"f a7rv -f a7rvl \ {p.p3){pi.p 2 ) , . 

* ^ - to4 + to 2 ,, - to2 J r< + (p - p 3 ) 2 + lP2 " P3 ' 

Macular = -iV^TWW r ^J** ^ + (<J <- / ) (B2) 
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The p contributions are: 



M 



2g*F^(m*-m*) 



m 2 + (p-p 3 y 



Aa-ycosQpgpw P2 ■ (gg - Pi) 



(p - P2) 



+ (P2 -^P3) 



(B3) 



The two point vertices are: 



C7T 

7TT7 



TTTj 



S~17T 
TTTj f 



Aacr 



1 cos9 p S' 6Ay A' 2 cos# p 4y a p cos6> p 



F 2 

Trp 


F 2 + 

Trp 


q 2 F 2 

•y Trp 


8y cos0 p 5' 


64yA' 2 cos6» p 


4ya p cos^p 


F 2 

Trp 


F 2 + 

Trp 


9 2 F 2 p 


8y sin8p 5' 
F 2 

Trp 


My A' 2 sin(9 p 
F 2 

■np 


Ay a p sm9p m 2 
9 2 F 2 p 


Sy s'mOp 5' 


QAyX' 2 sm9 p 


Ay otpSmOp m 2 , 



Tip 



F 2 



9 2 F 2 



Tip 



2y (b + d) sm0 s — V2y d cos9 s 
—2y(b + d) cos9 s — \/2y d sin^ s 



(B4) 



Notice that the superscript on C indicates which of the two particles involved in the AI = 1 transition is on-shcll; 
this only affects the a p term which has derivative coupling. 

It is not difficult to verify that the 7/ — > 7r°7r + 7r~ amplitude may be gotten from the one above by simply making 
the interchanges 



rj *-* r) , cosOp <-> sin# p 



(B5) 



everywhere in Eqs. (B1)-(B3). This should not be done in Eqs. (B4) since changing, for example, C% v to , 
accomplishes the desired result automatically. 
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